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ABSTRACT 

We report the late-time evolution of Type Hb Supernova (SN Hb) 2013df. SN 20I3df showed a 
dramatic change in its spectral features at ^ 1 year after the explosion. Early on it showed typical 
characteristics shared by SNe Hb/Ib/Ic dominated by metal emission lines, while later on it was 
dominated by broad and flat-topped Ha and He I emissions. The late-time spectra are strikingly 
similar to SN Hb 1993J, which is the only previous example clearly showing the same transition. 
This late-time evolution is fully explained by a change in the energy input from the ^®Co decay to 
the interaction between the SN ejecta and dense circumstellar matter (GSM). The mass loss rate is 
derived to be ^ (5.4 ± 3.2) x yr“® (for the wind velocity of ^ 20 km s“®), similar to SN 

1993J but larger than SN Hb 201 Idh by an order of magnitude. The striking similarity between SNe 
2013df and 1993J in the (candidate) progenitors and the GSM environments, and the contrast in these 
natures to SN 2011dh, infer that there is a link between the natures of the progenitor and the mass 
loss: SNe Hb with a more extended progenitor have experienced a much stronger mass loss in the final 
centuries toward the explosion. It might indicate that SNe Hb from a more extended progenitor are 
the explosions during a strong binary interaction phase, while those from a less extended progenitor 
have a delay between the strong binary interaction and the explosion. 

Subject headings: Circumstellar matter - stars: mass-loss - supernovae: individual: SNe 1993J, 
201Idh,2013df 
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1. INTRODUCTION 

Evolution of a massive star in the final stage toward 
the supernova (SN) explosion is one of the main nn- 
resolved issues in modern stellar astrophysics. Evolu¬ 
tionary paths to SNe Ilb/Ib/Ic, called stripped-envelope 
(SE-) SNe, have been actively debated - it is widely ac¬ 
cepted that the progenitors of SE-SNe have lost most 
SNe Ilb) or near ly all (SNe Ib/c) of an H-rich envelope 
Filippen^ 119971) . and the unresolved mass loss mech¬ 
anism should be a direct consequence of the progeni¬ 
tor evolution (iNomoto et~anil99t^ IWooslev et al.l[l994 
iGeorgy et al.ll2012l: iBenvenuto et al.ll2013[) . SNe lib are 
in particular important objects in this respect. They 
have a little but non-negligible amount of H-rich en¬ 
velope (< IMq) present at the time of the explosion 
(|Shigevama et al.l 11994 iBersten et al.l l2012i) , represent¬ 
ing a link between SNe II (with a red supergiant progen¬ 
itor) and SNe Ib/c (with a Wolf-Rayet like progenitor). 
Direct progenitor search in deep pre-explosion images 
points to a diversity in the progenitor, highlighted by dif¬ 
ferent extension of the pro genitor envel ope for SNe lib 
1993J (jMaund et al.ll2004 ) and 2011dh (|Van Dvk et al.l 
l2013f) : while both were nnprecedentedly identified as Yel¬ 
low Supergiant (YSG) stars, their sizes were significantly 
different with R ^ 6 OOi ?0 for the former and ~ 200i?o 
for the latter. 

SN lib 2013df was discove red on June 7.9. 2013 
UT soon after the explosion (jGhiabattari et aP 120131: 
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ICenko et alJ l2013fl i n NGC4414 at ~ 16.6 Mpc 

( Freedman et al.ll200l[) . A progenitor candidate has been 
identified in pre-SN HST images, whose nature is quite 
similar to the progenitor of SN 1993J, interpreted as 
a gia nt with a radius of 545 ± 65i?0 (jVan Dvk et al.l 
l2014ll . Pr operties as an SN are also similar to those of 
SN 1 9 93J (iMorales-Garoff o lo et al.ll20l4 l\^n Dvk et al.l 
I20l4 iBen-Ami et al.l l^isfh both showed similarly a 
bright ‘early-cooling’ phase, and similar ‘peak’ light 
curves (LGs) and spectral evolutions around the maxi¬ 
mum light. Throughout this paper, we refer the emis¬ 
sion from the adiabatically expanding SN ejecta follow¬ 
ing the shock breakout as the ‘early-cooling emission’; 
the strength and duration of t his emission are sen sitive 
to the progenitor radius (e.g., iNakar fc Piroll20f^ . and 
it lasts for less than ^10 days after the shock breakout 
for SNe Ilb. The phase is followed by the emission pow¬ 
ered by radioactive decays of ®®Ni to ®®Go and then to 
^®Fe, corresponding to the main part of the LG which is 
peaked at ~ 20 days after the breakout for SNe lib. This 
‘peak’ LG and spectral properties are mostly determined 
by the properties of the SN ejecta rather than those of 
the progenitor. While the ejecta mass has been estimated 
slight ly smaller in SN 2013df (jMorales-Garoffolo et al.l 
l2014f) . these two SNe seem to be almost a twin in terms 
of the progenitor and the SN properties. A contrast 
of these SNe to another well-studied SN lib 20IIdh is 
interesting - while showing the similar peak LGs and 
spectral properties (thus si milar SN properties), its pro¬ 
genitor was less extended ((Van Dvk et al.l[20I3ll as was 
consistent with the faint early-phase cooling emission 
(|Bersten et al.|[20T^ . 

In this paper, we report the late-time optical evolu¬ 
tion of SN 2013df up to > 600 days after the explosion. 
§2 presents observations and data reduction. §3 shows 
that the spectral and light curve evolutions are found 
to be dramatically changed in the late phase, showing 
clear signatures that the SN is now powered by the in¬ 
teraction between the SN ejecta and circumstellar matter 
(GSM). Among SE-SNe, such a transition was observed 
previously only for SN 1993J, and SN 2013df is the sec¬ 
ond example. Indeed, the present study reveals that SN 
2013df is a twin of SN 1993J, not only in the early phase 
but also in the late phase. This again provides a striking 
difference to SN 201 Idh, for which a strong interaction 
was not observed in the late-phase. In §4, we analyze the 
late-time spectra and light curves, deriving properties of 
GSM, and discuss a relation between the natures of a 
progenitor and the mass loss. The paper is closed in §5. 

2. OBSERVATIONS AND DATA REDUCTION 
2 .1. Spectra 

Our spectroscopic observations were performed on 
MJD 56648.6 (22 December 2013; day 199) and on MJD 
57075.4 (22 February 2015; day 626) with the 8.2m 
Subaru telescope equipped with the Faint Object Gam - 
era and Spectrograph (FOGAS; iKashikawa et al.l[2002[l . 
Additional two spectra were obtained on MJD 56604.5 
(8 November 2013; day 155) and on MJD 56783.3 (6 
May 2014; day 333) with the 6.5m MMT tele scope using 
the Blue Channel instrument (jSchmidt et al.l[T989tl . The 
dates indicated here are measur ed from the well-defined 
explosion date (MJD 56449.5: [Morales- Garoffolo et al.l 


I20l4 IVan Dvk et al.l 120141 , and this definition is used 
throughout the paper. 

The setup for the Subaru/FOCAS observations is the 
following: We used 0.8” slit and the B300 and R300 
grisms on day 199 covering 3,650-10,OOOA, and 0.8” off¬ 
set slit and the B300 grism on day 626 covering 4,700- 
9,000A. The slit direction is set at PA = 0 using an 
Atmospheric Dispersion Corrector (ADC). The spectral 
resolution is ~ 500. The exposure times are 600 s (in 
AB) and 3600 s (in ABBA), respectively. For flux cali- 
brati on and telluric absorption correction, Feige 34 (lOkd 
1199011 was observed in both nights. For MMT Blue 
Channel observations we used a 1” slit and the 300 line 
mm“^ grating covering 3,200-8,500A. Both observations 
were at the parallactic angle. The spectral resolution is 
R = 740. The exposure times are 600 s (single) and 
2 X 1200 s, respectively. For flux calibration and t elluric 
abs orption correction, s t andar d stars from lOkd (|1990ll 
and iMassev fc Gronwalll (|1990[ 1 were observed. 

The spectra are reduced following standard proce¬ 
dures with IRAF/PyRAF,^^ including bias subtrac- 
tion, flat fielding, c osmic ray rejection (with LAcosmic; 
IVan DokkumI l2001h . sky subtraction, ID spectral ex¬ 
traction, wavelength calibration using ThAr or HeNeAr 
lamps, and flux calibration. The wavelength solution 
was checked with the atmospheric sky emission lines, and 
when necessary a small shift was further applied. 

2 .2. Imaging and Photometry 

For photometry, we obtained V (20 s) and R (20 s) 
images on day 199 and V (120 s) and R (180 s) im¬ 
ages on day 626 with Subaru/FOC AS. On day 19 9, the 
standard stars around PG1525-071 ()Landoltlll992ll were 
observed. We also obtained the SDSS r-band images at 
three epochs (300 s for each), on day 262, 298, and 348, 
using MMTCam. 

PSF photometry was performed on day 199 using 
the images around the PG1525-071 under photometric 
condition. This results in V = 19.2 ±0.1 mag and 
R — 18.5 ± 0.1 mag. The photometry on the MMTCam 
images results in r = 19.75 ±0.10 (day 262), 20.08 ±0.08 
(day 298), and 20.59 ± 0.09 (day 348). The MMTCam 
images were calibrated using observations of SDSS field 
stars. Since the SN was found to be still relatively bright 
in these epochs, we have not subtracted a reference image 
for photometry with these images. 

The SN has faded substantially on day 626, but it was 
clearly detected in both V and R bands (Fig. 1). Since 
the SN magnitude is now comparable to the background 
which shows substantial fluctuation around the SN site, 
the photometry on day 626 requires a careful analysis. 
First we performed relative photometry between the im¬ 
ages on day 199 and 626, then we subtracted the pre-HST 
WFPC2/F555W image on 29 April 1999 (Obs ID: 8400, 
PI: Keith Noll; see also IVan Dvk et ^l2014f) from our 
P-band image on day 626. Thereby we have obtained 
V ^ 22.6 mag. This was used to anchor the flux scale in 
the spectrum, then the R and /-band magnitudes were 

IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. PyRAF is a product of the Space 
Telescope Science Institute, which is operated by AURA for NASA. 
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derived as i? ~ 22.0 mag and I ~ 22.5 mag from spec¬ 
troscopic photometry. Details on the photometry will be 
presented in a separate paper (Maeda et ah, in prep.). In 
any case, the uncertainty of the flux level even by a fac¬ 
tor of ^ 2 would not affect our analyses and conclusions 
in this paper. 


3. RESULTS 

Figure 2 shows the time sequence of late-time spec¬ 
tra of SN 2013df from day 155 to day 626. Fig¬ 
ure 3a shows the spectra of SN 2013df at day 199 
and day 626, as compa red to those of SNe Ilb 1993J 
(iMatheson et al.l l2000al lbll and 201 Idh (|Shivvers et al.l 
1201 .311 at similar epochs. The comparison spectra are 
downloaded fro m The Weizmann inter active supernova 
dat a repository dOfer fc Gal- Yarn! 201211 .^° As discussed 
by iMorales-Garoffolo et al.l (12014 1. the spectra of SN 
2013df at < 200 days since the explosion are overall typ¬ 
ical for SNe Ilb/Ib/Ic, dominated mostly by forbidden 
lines of metals and some allowed transitions. There are 
two especially interesting features: (1) The ratio of [01] 
6300, 6363 to [Ca II] 7291, 7324 is low, and (2) there is an 
excessive emission on the red side of the [O I]. The feature 
(1) can be interpreted as a relatively low-mass prog eni- 
tor (jFransson fc Chevaliedll98^lMaeda et al.l[^07al1 . in¬ 
dicating the importance of the binary interaction to get 
rid of most of the hydrogen envelope before the explosion 
to become an SN lib. The feature (2) i s quite common in 
SNe lib, as exemplified by SNe 1993J IM^ 

Matheson et al.ll2000al fOl and 201 Idh (|Sahu et al.ll2013l 


Shivvers et al. 'l2013ll. The similar feature is also seen 


in soine SNe Ib/c (ISoller man et al.l Il998l iMaeda et al.l 
l2007bl : iTaubenberger et al.l 20111 : iFolatelli et ahll^Ma l. 
A straightforward identification is Ha, but the solid iden¬ 
tification is difficult - the featu re could be attributed t o 
[N II] 6548,6583 at ~ 200 days (j.lerkstrand et al.ll2015h . 

A dramatic change was observed in the spectrum taken 
on day 626. There is little trace of metal emission lines 
anymore, but the spectrum is characterized by broad 
emission features including those around ^ 6, 550A and 
^ 5, 800A. Indeed, a resemblance of this spectrum to 
that of SN 1993J at a similar epoch is very striking (Fig¬ 
ure 3b). At such a late epoch, a canonical Co-powered 
nebular emissio n model is unable to ex plain these broad 
emission lines (jJerkstrand et al]l2015ll . and there is al¬ 
most no doubt that these features are Ha and He I 
5876/Na ID (plus [N H] on the blue side), respectively. 
This spectrum shows that the interaction between the SN 
ejecta and CSM is now a predominant power source. The 
spectrum on day 333 nicely shows the development of 
Ha and the onset of a transition from the ®®Co-powered 
emission to the SN-CSM interaction-powered emission 
(Figure 2). This is t he second example, wi th SN 1993J 
as the first example ([Matheson et al.lT2000al lbh . where an 
SN Hb/Ib/Ic shows such a transition. A marginal detec¬ 
tion of a broad Ha was also reporte d for SN 201 Idh at a 
similar epoch ([Shivvers et al.l [201^ . but the strength of 
the Ha does not compare to those found in SNe 2013df 
and 1993J (Figure 3a). 

Figure 4a compares the profiles of Ha emission seen 
in SNe 2013df and 1993J at > 600 days. First, for both 
SNe 2013df and 1993J, the Ha shows a boxy profile with 
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a steep cut-off on both sides. This indicates that the 
main part of the emission is originated in a thin shell 
with the velocity of ^ 10,000 km s“^ (see the model 
curves in Figure 4a). There is however noticeable dif¬ 
ference in details, where SN 2013df shows an excessive 
emission near the rest wavelength but blueward, while 
SN 1993J shows a suppression in the same wavelength 
region. A similar behavior is marginally visible in both 
blue and red edges, where SN 2013 df shows an excess (see 
also lMorales-Garoffolo et al.ll2014fl . The difference in the 
line profiles is probably caused by a different structure 
and/or viewing direction (see below for further support). 
If it is put into a unified context, a density enhance¬ 
ment in the equatorial direction, with SN 2013df viewed 
from the polar direction and SN 1993J from the equa¬ 
torial direction, could explain the main difference seen 
around the rest wavelength. There are another couple of 
possible, but less likely, explanations for this relatively 
narrow, central component in SN 2013df: (1) This might 
reflect an unshocked GSM ionized by the high-energy ra¬ 
diation from the shocked region. The velocity is how¬ 
ever too high (^ 2,000 km s“^) . (2) This r night reflect 
swep t-up materials from a companion star ([Hirai et al.l 
I 20 TI. analogous to what has been searched for in SNe 
la (iLundavist et ^120131 : [Maeda et al.ll20I4Hl . If we ap¬ 
ply the model prediction for SNe la, the feature in SN 
2013df seems too strong to be interpreted this way. 

Figure 4b shows the feature at 5, 8 OOA. The red part 
of the feature matches well to Ha assuming this feature is 
either He I 5876 or Na ID. This supports an idea that this 
feature is also powered by the SN-GSM interaction. The 
blue part shows some excessive emissions, which might 
be contaminated by [N H] and [O I], either from the 
SN ejecta, the interacting region, or the unshocked but 
irradiated CSM. 

Figures 4c shows the evolution of SN 2013df from 
day 155 to day 626. If we scale the spectrum at day 
626 multiplied by a factor of ^ 20, 8, and 2.5, the 
emission at the red side of the [O I] as seen in the 
spectra at day 155, 199, and 333 matches to the Ha 
profile very well. This suggests that this feature seen 
at as early as day 155 was already dominated by Ha . 
This was speculated by IMorales-Garoffolo et al.l ([201411 , 
but the late-time spectrum is necessary for the conclu¬ 
sive i dentification, given po ssible alternative explana¬ 
tions ([Jerkstrand et al]l2015l1 . Since the profile of Ha 
does not show substantial evolution from day 155 to day 
626 (except for the high-velocity edge; see below and 
§4.1), the line profile should reflect a geometrical nature 
of the emitting region, rather than the contributions by 
other lines, e.g., [N H]. 

Figure 4d compares the red part of the Ha profile on 
day 199 and day 626. While the overall profile is kept 
nearly the same, there is a decrease in the velocity cor¬ 
responding to the wavelength at the cut-off, namely the 
velocity of the emitting shell. This is qualitatively in 
agreement with the expectation from the SN-CSM inter¬ 
action, where the shell created at the interacting region 
should be decelerated as a function of time. This will be 
further discussed in §4.1. 

Figure 5 shows long-term V and i?-band light curves 
of SN 2013df as compared to those of SNe 1993J and 
201 Idh. In this figure, we have converted the SDSS r- 
band magnitudes obtained by the MMTCam to the R- 
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band maKnitudes assumiiiK a cons tant offset of —0.15 
mag (jMorales-Garoffolo et al.l[2014l ^. Up to ^ 200 days, 
the three SNe show mutually similar evolutions, sug¬ 
gesting that the ejecta mass (Mej) and the kinetic en¬ 
ergy (Sk), or specifically the opti cal depth to y-ray s 
from ®®Co as scaled by M^-/Ek (iMaeda et al.l 120031 1 . 
are similar; in this sense, the ejecta mass obtained by 
iMorales-Garoffolo et al.l (I2014I1 . Mgj = 1.2 ± 0.2Mq and 
Ek = 0.8 ±0.4 X 10®^ erg for SN 2013df, as compared to 
Mej = 2.7±0.8Mf7, and.E K = 1.2 ± 0.2 x lO’^^ erg for SN 
1993J (lYoung et al.lll99^ . is probably underestimated. 
On day 626, SN 2013df turns out to be brighter than 
SN 2011dh by ^ 2 mag in the U-band and 2.5 mag 
in the i?-band (when normalized at their peak magni¬ 
tudes) , as is consistent with the power source dominated 
by the SN-CSM interaction. This suggests that the main 
energy source had changed at ^ 1 years after the explo¬ 
sion. Indeed, the i?-band light curves of SNe 2013df and 
1993J both deviate from that of 2011dh at ^ 200 — 300 
days (as is also seen in the U-band light curve of SN 
1993J), exceeding the magnitudes of SN 2011dh at the 
corresponding epochs. Together with the spectral sim¬ 
ilarity, this also suggests that the SN-CSM interaction 
property is very similar between SN 2013df and 1993J. 
We note that the similarity in the SN-CSM interaction is 
also supported by an independe nt analysis of radio and 
X-ray properties of SN 2013df (jKambie et al.l [2?)T^ as 
further mentioned in §4. 

Currently SN 2013df is much brighter in the V and R- 
bands than the expectation from the ®®Co-decay heating, 
and this would raise a hurdle for searching for disappear¬ 
ance of the progenitor candidate a nd/or for emergence 


Van Dvk et al.l 120131: 

Folatelli et al.l l2014bl: iFox et al.l 

20141). IVan Dvk et al. 

(|2014l) identified the candidate 


progenitor as a giant with i?eff = 545 ± 65i?Q with 
V = 24.5 mag and / = 23.1 mag {R ~ 24 mag if simply 
interpolated). If there would have been only the ®®Co 
decay power (following the light curve of SN 20IIdh), 
the V and i?-band magnitudes should have been already 
comparable to or even fainter than those of the progeni¬ 
tor candidates on day 626. Due to the additional power 
input by the interaction, the epoch when the SN becomes 
fainter than the progenitor candidate should be delayed. 

In Figure 5b, we show the contribution from the in¬ 
teraction power in the i?-band light curves of SN 2013df 
(denoted by star symbols) estimated as follow: Assuming 
that the Ha luminosity follows the interaction power, we 
estimated the interaction power at the epochs where the 
spectra are available (Figure 4c). Then we assumed that 
the i?-band magnitude at day 626 is dominated by the 
interaction, and normalized the interaction power at this 
epoch. The same evolution is assumed for the U-band. 
Avoiding possible overestimate of the interaction power 
in the relatively early phase (i.e., day 155, through the 
contribution by [N II] and continuum), we fit the inter¬ 
action power thus estimated on days 199, 333, and 626. 
With only three points the fit is uncertain, but we should 
be able to provide some meaningful estimate on when the 
confirmation of the progenitor is possible. If we assume 
an exponential decline in luminosity then the fit results 
in L oc while if we assume a power law then 

L oc In the exponential case, the SN will become 


fainter than the progenitor candidate at ^ 1,000 days af¬ 
ter the explosion both in V and i?-bands, i.e., February - 
March 2016. In the case of the power law decline, it will 
further be delayed to ~ 1,600 days after the explosion, 
i.e., October - November 2017. We suggest continues 
follow-up observations of SN 2013df to refine the predic¬ 
tion. This point will be further discussed in a separate 
paper (Maeda et ah, in prep.). 

4. DISCUSSION 

4.1. Analyses of the Ha and CSM around SN 2013df 

The observed broad and flat-topped Ha emission is 
naturally expected as arising from the SN-CSM inter¬ 
action. The interaction creates a reverse shock (RS), 
contact discontin uity (CD), and forward shock (FS) 
(jChevaliei]ll982all . The Ha emission can be originated 
in the unshocked SN ejecta irradiated by X-ray photons 
created at the RS: Ionizations induced by high-energy ra¬ 
diations are followed by recombinations. This is also an 
interpr etation for Ha emission o bserved in middle-aged 
SNe H (iMilisavlievic et ani2012ll . Alternatively, it could 
be contributed by the RS region if the temperature in 
the RS region is decreased to ^ 10,000 K. The latter 
situation is less likely as the RS is typically considered 
to have much hi gher temperature to emit X-rays for SNe 
Hb/Ib/Ic (e.g., iMaeda et al.ll2014ah . and we note that 
the FS region would have even much higher temperature. 
Indeed, analysis of the X-ray emission from SN 2013df 
does indicate the temperatures in the RS and FS regions 
as > lO^K and > 5 x 10®K, respectively (jKambie et al.l 
1201511 . In either case of the unshocked ejecta or the RS 
region, the outer edge of the Ha-emitting region can be 
identified as the CD. The emitting region may well be 
identified as nearly a whole H-rich envelope in the un¬ 
shocked ejecta scenario while as a part of the region in 
the RS scenario. 

For the distance of 16.6±0.4 Mpc, the Ha luminosity is 
~ 8.8x10^® erg s“^ on day 199 and ~ 1.1x10®® erg s“® on 
day 626, with E{B — V)= 0.1 mag and Ry = 3.1. This 
indicates that the Ha luminosity is similar to, or a bit 
stronger than, SN 1993J, for w hich it was derived to be ^ 
1.26 X10®® erg s“® on day 367 (iPatat et al.lll995ll . At day 
626, the velocity of the shell is ~ 10,000 km s“® (Figures 
4a and 4d). An uncertainty of ^ 5% in the shell velocity 
is mainly from the uncertainty in defining the continuum. 
From a comparison between the Ha profiles between day 
199 and day 626 (Fig. 4d), we observe the velocity has 
decreased by ^ 5 — 10% (where the main uncertainty 
is from the estimate of the shell velocity as mentioned 
above). The emitting region should be confined within a 
thin shell with ARsh/Rsh ^ 0.2 as constrained by the line 
profile (Figure 4a). Within the fitting error, Ai?sh/Rsh 
0.1 is acceptable, but ARgh/Rsh ^ 0.3 is rejected. 

As for the ‘He F feature, the flux in the He I/Na ID at 
day 626 is by a factor of ^ 3.2 smaller than Ha (excluding 
the excessive emission blueward: Figure 4b). Attributing 
this feature to He I 5876, we can derive the abundance 
in the emitting region. Assuming the Case B recombi¬ 
nation both for Ha and He I, the relative flux will be 
LifAa)/LilAel) ^ 1.3n(H+)/n(He®’) for a wide range of 
tempe rature (5,000 - 20 OOOK) where n’s are number den¬ 
sities (|Osterbrocklll989| l . For the observed ratio of 3.2, it 
requires n(H+)/n(He~'") ^ 2.5. Assuming that all the hy- 
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drogen and helium are singly ionized, the emitting region 
must be composed of 70% of H and ~ 30% of He in 
number, or ^ 40% of H and ^ 60% of He in mass. Typ¬ 
ical models for SN Hb adopt the He-rich H/He envelope 
with the mass fraction of He larger than that o f H (e.g., 
V ^ 0.8 for SN 1993J: IShigevama et al.lll994li . and our 
derived abundance in the emitting region is consistent 
w ith the expe c tation . 

iPatat etahl (|1995ll constrained the mass of the enve¬ 
lope hydrogen for SN 1993J on day 367. Since the derived 
geometry of the emitting region for SN 2013df is simi¬ 
lar to the situation they considered, we use their equa¬ 
tion 2 to estimate the envelope hydrogen mass for SN 
2013df. In doing this, we have an additional constraint 
for the electron density as 1.4n^ (see above). Adopt¬ 
ing the velocity at the CD as t4h = 10,000 km s“^ and 
the luminosity at day 626, we derive M(H) ^ O.25M0 
at day 626 for the temperature of 10,000K. Adopting 
14h = 11,000 km s“^ (see above; Figure 4d), we de¬ 
rive M(H) ~ O.14M0 at day 199. We note these two 
estimates are consistent within the uncertainty in the 
temperature of the emitting region: If the temperature 
on day 199 would be ~ 20, OOOK (i.e., the interaction 
power is higher in the earlier phase), then the estimate 
goes up to M(H) ~ O.28M0. Therefore, we conclude 
that the mass of emitting hydrogen is M{H) ^ 0.2Mq. 
The lack of the evolution favors the scenario where the 
emitting region is the unshocked ejecta: If the Ha emit¬ 
ting region is the RS, then the mass should have been 
increased by a factor of more than three from day 199 
to 626. The mass we derived is therefore likely repre¬ 
sentative of the whole H-rich envelope. If we assume the 
composition we derived for SN 2013df is also a pplicable 
to SN 1993J, the estimate bv iPatat et al.l (jl995l) becomes 
M(H) ~ O.14M0. Therefore, the mass of the envelope 
hydrogen is similar between SN 2013df and 1993J, with 
SN 2013df probably larger by a factor of ^ 1.4. By com¬ 
bining H and He, we derive the total mass of the emitting 
region as ~ 0.5Mq. 

From the profile of Ha emission, we estimate that 
the expansion velocity of the shell decreases follow¬ 
ing t4h oc where /3 ~ — 0.04 to ~ (Figure 

4d). The self-similar solution (|Chevalieijll982ali predicts 
P = {3 — s)/{m — s), where m and s are power law in¬ 
dexes in the density distribution as a function of veloc¬ 
ity or radius for the SN ejecta and the CSM, respec¬ 
tively. Adopting s = 2 as expected for the steady-state 
mass loss which is (roughly ) applicable t o SNe 1993J 
(iFransson fc Chevalied fl9^ and 2011dh (jMaeda et al.l 
l2014al) . we derive m ~ 15 — 27 which takes into ac¬ 
count the uncertainty in measuring the deceleration of 
the shell velocity. While measuring the exact value of 
m is difficult due to the sensitive dependence of m to 
P, the outer ejecta structure of SN 2013df seems to be 
steep. This steep density gradient is difficult to reconcile 
with a stand ard model of the SN dynamic s which predict 
TO ^ 7 — 12 (IChevalier fc Franssonl [200^ , but is similar 
to wh at have been derived for SNe 1993J (IFransson et al.l 
Il996l) and 201 Idh (|Maeda et al.l 1201^ . This indicates 
that our understanding of the outermost density struc¬ 
ture in the H-rich env elope of the SN Hb progenitors is 
still missing (see e.g., iNakar fc Piroll2014h . and the fact 
that the steep density gradient is generally found in the 


well-studied SNe Hb should provide a hint in understand¬ 
ing the still unresolved final evolution of the progenitors 
toward SNe Hb. 

Applying the self-similar solution of lChevalied (|1982ali 
for the velocity a t the CD, specifically using eq. 5 of 
iPatat et"^ (jl995l l. we can derive the CSM density, i.e., 
A* as defined by pcsM = 5 x 10^^A*r“^ (i.e., A^, = 100 
for the mass loss rate of ~ 10~^Mq yr“^ and the wind 
velocity of ^ 10 km s~^). For a range of to ~ 15 — 
27, Mej = 2Mq (iMorales-Garoffolo et al.ll2014f) . and the 
characteristic ejecta velocity of 10,000 km s“^, we derive 
A* ^ 110 (for TO = 27) and A* ^ 430 (for to = 15) 
for SN 2013df. Therefore, we conclude that the CSM 
around SN 2013df is dense with A* ^ 110 — 430. If this 
CSM density is converted to the mass loss rate under the 
assumption of a steady-state mass loss (which roughly 
applies to SNe 1993J and 2011dh: iMaeda et al.ll2014^ 
and the wind velocity of ~ 20 km s“^, the mass loss 
rate in the final ~ 800 years before the SN explosion 
is (5.4 ± 3.2) X 10~^Mq yr“^. Throughout the paper, 
we adopt 20 km s“^ as our fiducial value for the SN Hb 
progenitor wind velocity; this is at the low end as inferred 
for YSGs (20 — 100 km s“^; ISmithl[2?)l^ . and the wind 
velocity of ~ 10 — 20 km s“^ have be en assumed for 
the p r ogenitors of SNe 199 3J and 201 Idh (IFransson et al.l 
119961: iMaeda et al.ll20ll^ . We note that in any case 
the uncertainty in the wind velocity by a factor of a few 
would not affect the conclusions in this paper. Moreover, 
the derived mass loss rate is simply scaled linearly with 
the assumed wind velocity. 

The derived CSM density is also largely consistent with 
X-ray properties of SN 2013df. SN 2013df was brighter 
than SN 2011dh in X-rays by a fa ctor of ~ 10 (|Li fc Kond 
l2013l) . Figure 9 of IMaeda et al.l (l2014al l. taking into ac¬ 
count of the cooling and absorption in the RS region, 
shows that the predicted X-ray luminosity is roughly con¬ 
sistent with that of SN 2013df for A* ^ 200—400. Indeed, 
from a detailed an alysis of the ra d io an d X-ray proper¬ 
ties of SN 2013df, iKamble et al.l (l2015li independently 
derived the CSM density which is roughly in agreement 
with the estimate mentioned here. 

The mass in the RS region (i.e., the mass in the swept- 
up ejec ta) can be estima ted for given A* and to (and 
s = 2). IChevalieil (jl982bl l showed the scaling relation as 
M^is/Mps = (to — 4)/2, where Mrs and Mrs are the 
masses in the RS and FS regions. This formula overes¬ 
timates the value found in a more detailed self-similar 
analysis of IChevalieil (jl982al) by ^ 30%. Using these 
relations, we estimate Mrs 0.28Mq for to = 15 and 
Mrs ~ O.15M0 for TO = 27. Note that Mrs is not 
sensitive to the exact value of to. Taking into account 
this uncertainty, Mrs does not seem to be sufficient to 
explain the mass of the Ha-emitting region (~ 0.5Mq), 
adding another support for the interpretation that Ha is 
produced by recombinations within the unshocked ejecta. 
The temperature behind the RS is probably too high to 
contrib ute to the Ha (e.g. , 1 — 3 keV derived for SN 
2011dh: IMaeda et al.ll2014ah . and then we derive the en¬ 
velope mass as ~ 0.7Mq by adding this mass to that 
of the Ha emitting region as derived above. Note that 
this estimate is fully independent from LC models us¬ 
ing the early-phase LC at the fir st few days which i s 
yet to be fully developed (e.g., iNakar fc Pirol [201^ . 
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and our late-time observation provides a good calibra¬ 
tion for such models. Our result indicates similar struc¬ 
tures in the progenitor envelopes between SNe 2013df 
and 1993J - the H-rich envelope mass was estimated to 
be 0.3 — IMq through the LC models for SN 1993J. 
Note that the H-rich envelope masses in SNe 2013df and 
1993J are substantially larger than that o f SN 2011dh, 
which was constrained to be ^ O.IM^ (lArcavi et al.l 
I2OIII: iBersten et al.l 120121 : iMaeda et al.l boidall . This 
again points to the difference between SN 2013df (and 
SN 1993J) and SN 2011dh. 


4.2. A Link between the Natures of Progenitors and 
Their Mass Losses 

A sample of well-studied SNe Ilb to date, including the 
pre-SN progenitor (candidate) detection, prompt multi¬ 
wavelength follow up in a few days since the explosion, 
and a long term follow-up covering more than 500 days, 
composes of only three SNe lib (1993J, 2011dh, and 
2013df). However, the similarity and difference among 
the sample are striking. It has been shown that these 
three are similar in the SN ejecta properties as inferred by 
their maximum-light prope rties (jMorales-Garoffolo et al.l 
I20l4 IVan Dvk et al.l 120141 !. The nature of the progen¬ 
itor seems to be divided into two classes, a more ex¬ 
tended progenitor (SNe 1993J and 2013df) and a less 
extended one (SN 201 Idh), as inferred bot h by the pro¬ 
genitor (candidate) in the pre-SN images (iMaund et al.l 
|2004 I Van Dvk et al.l l2013l I20l4 and by the behavior 
in th e LC in the first few d ays ( Morales-Caroffolo et al.l 
I20l4 IVan Dvk et al.ll20i4l . This classification and di¬ 
versity are further extende d by adding UV pro perties 
to the list of observables (|Ben-Ami et al.l l2015l) . We 
find that this classification further applies to the late¬ 
time behavior at > 600 days since the explosion. While 
the sample is still limited, the similarities between SNe 
1993J and 2013df, as well as the contrast to SN 2011dh, 
are so striking, and we speculate that these two classes 
might reflect distinctly different paths in the evolution 
toward an SN. This might be related to the ‘extended’ 
and ‘compact’ classification fo r SNe Hb as proposed by 
IChevalier fc Soderberd (l2010h m ostly from their radio 
properties fsee also lKamble et al.ll2015l) . 

In Figure 6, we show comparison of SNe Hb in prop¬ 
erties of the early-phase ‘cooling’ emission / progenitor 
radius and the CSM density / the mass loss rate in the fi¬ 
nal centuries toward the SN explosion. The properties in¬ 
ferred for SNe Ib/c are also indicated, while for SNe Ib/c 
the derivation of the CSM density and thus the mass loss 
rates is no t as robust as that for the well-studied SNe Hb 
(see, e.g., iMaedal[20121 for the uncertainty). In convert¬ 
ing the CSM densities to the mass loss rates, we assume 
the mass loss wind velocity of 20 km s“^ for SNe Hb and 
1,000 km s“^ for SNe Ib/c. 

There is a correlation between the strength of the 
early-cooling luminosity and the CSM density as shown 
in Figure 6a. SNe Hb 201 Idh and 2008ax have the 
CSM density at the high-end assumed for SNe Ib/c, 
both of which have bee n classified as ‘compact’ SNe Hb 
([Soderberg et al.l 12012!) .^^ SNe Hb 2013df and 1993J 
clearly have a larger CSM density than these SNe Hb 


Given the confirmed YSG progenitor of SN 2011dh, more 
appropriate terminology would be ‘less extended’. 


by an order of magnitude, and than typical SNe Ib/c 
by two orders of magnitudes. SN 1993J reached a peak 
in the early-cooling emission, as seen in the pseudo- 
bolometric optical-NIR LC, at ~ 3 — 4 days after the 
explosion. This is due to the combination of decreasing 
bolometric luminos ity and decreasing temperature (e.g., 
iNakar fc Pircill2M4) . As the early-cooling luminosity, we 
adopt this luminosity for SN 1993J, and the pseudo- 
bolometric optical-NIR luminosity (or upper limit) at 
the same epoch for the other SNe Hb (see, e.g., figure 
6 of iMorales-Garoffolo et al.ll20l4l . SNe Ib/c and ‘com¬ 
pact’ SNe Hb do not show signatures of the early-cooling 
emission (except for SN 2011dh thanks to a deep and 
prompt follow-up), but deep limits have been obtained 
for some SNe Ib/c (e.g., ICao et aHl2013D . No progeni¬ 
tor candidates have been reported for SNe I b/c except 
for a possible WR progenitor of iPTF13bvn (ICao et al.l 
1201311 .^^ In any case, a giant progenitor has been rejected 
for some SNe Ib/c (e.g., see lSmarttll2009l : lEldridge et al.l 
IMl, for a review). As expected, the relation is also 
clearly seen if we use the progenitor radius instead of the 
early-phase cooling luminosity (Figure 6b). This find¬ 
ing is consistent with the correlation between the prop¬ 
erties of early-cooling phase and UV properties around 
maximum-light if the lat ter is interpreted as a trace of 
the SN-CSM interaction (|Ben-Ami et al.l[2015lj . 

While there seems to be a well-defined relation between 
the progenitor radius and the CSM density, it becomes 
less obvious if the CSM density is converted to the mass 
loss rate assuming typical wind velocities for a WR star 
(^ 1,000 km s“^) for SNe Ib/c and a YSG star (~ 20 
km s“^) for SNe Hb (see §4.1), as shown in Figure 6c. 
We show this for the purpose of demonstration since the 
mass loss rate is more directly linked to the progenitor 
evolution than the CSM density, while the CSM density 
is more directly linked to the observations; the inferred 
mass loss rates here suffer from the uncertainty by a fac¬ 
tor of a few in the wind velocity. In any case, our argu¬ 
ment below is based on the difference in the mass loss 
rate by an order of magnitude, and thus the uncertainty 
in the wind velocity would not affect our conclusions. 
For SNe Ib/c this is derived for the radio observation 
at < 100 days after the explosion corresponding to the 
final < 10 years of the mass loss toward the explosion, 
and for SNe Hb the derivation is based on the observa¬ 
tions at ~ 500 days corresponding to ^ 1000 years before 
the explosion. The range of the mass loss rate inferred 
for SNe Ib/c seems to overlap with the range seen in 
SNe Hb. This mostly comes from the high wind veloc¬ 
ity for a hypothesized WR progenitor of SNe Ib/c, and 
from so me SNe Ib/c that sh ow a relatively high CSM 
density (IWellons et al.ll2M^ . In any case, if we focus 
on the well-observed SNe Hb, there is a relation (within 
the small sample) between the progenitor radius and the 
mass loss rate, divided into SNe 2013df and 1993J on one 
side (more extended and a larger mass loss rate) and SN 
2011dh on the other (less extended and a smaller mass 
loss rate). 

In sum, we suggest that there is a relation between the 

H owever, we note that it has been claimed by lEldridge et ahl 
<1201511 that the detected progenitor candidate is more likely a 
binary system c onsisting of lower mass stars than claimed by 
IGao et II2OI3II . 
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nature of the progenitor (i.e., in the radius and mass of 
the H-rich envelope) and the mass loss rate in the final 
centuries toward the explosion, for a giant progenitor of 
SNe Ilb. To clarify if this is a generic feature among SNe 
Ilb requires a larger sample, but we see this does apply 
to the well-observed three SNe Ilb 1993J, 2011dh and 
2013df: Progenitors with a more extended and massive 
H-rich envelope have larger mass loss rates in the final 
centuries toward the explosion. This is at odds with a 
naive expectation: If a progenitor had a larger mass loss 
rate, it would lose a larger amount of the H-rich enve¬ 
lope, contrary to what we observe. One possibility is that 
the relation reflects a possible causality between the two 
quantities: With such a high CSM density for SNe 2013df 
and 1993J (H* ~ 200), the electron scattering opacity in 
the wind region before the explosion, if the CSM was 
fully ionized (while unlikely for a giant progenitor), is 
estimated to be an order of unity near the vicinity of the 
surface of a progenitor. If this is the case, the progenitor 
radius might have been falsely identified in the photo¬ 
sphere within the wind region for SNe 2013df and 1993J. 
However, this does not explain another observational in¬ 
dication that we also see the difference between the two 
classes in the mass of the H-rich envelope. While further 
investigation is necessary, we regard this possibility as 
unlikely. 

Therefore, we suggest that this unexpected relation 
between the progenitor radius and the mass loss rate 
reflects a real difference in the evolutionary paths to 
these ‘two classes’ of SNe Hb. Given that the binary 
inter action is proposed to be a key for these SNe Hb 
('e.g.. lBenvenuto et al.ll2013ll and that the mass loss rate 
is affected s ubstantially by the binary interaction (e.g., 
iSmithI 1^014 for a review), this suggests that there is a 
difference in the binary interaction history for these two 
classes, reflecting differences in the evolution of the bi¬ 
nary separation and mass transfer history. The mass 
loss rate just before the explosion was found to be in- 
sufScient to produce a YSG progenitor for SN 2011dh 
(|Maeda et al.l l2014af) and a strong binary interaction 
phase well before the explosion (at least > 1,000 years 
before the explosion) is necessary. On the other hand, the 
derived mass loss rates for SNe 2013df and 1993J are in¬ 
deed enough to expel nearly all the H-rich envelope in the 
giant phase (if this mass loss rate would keep that high 
during a whole giant phase), while the binary interaction 
is probably required to drive such a large mass loss rate 
in the giant phase. Therefore, it seems that the timing of 
the binary interaction might be different between the two 
classes. As an emerging picture, we propose that SNe Hb 
with a more extended progenitor are explosions during a 
strong binary interaction phase (including any effects by 
a companion star such as irradiation), while those with a 
less extended progenitor have a time delay between the 
strong binary interaction phase and the explosion. A 
further investigation of a possible companion star for SN 
2013df, whose nature should be related to the binary in¬ 
teraction history, will be an important step to solve this 
issue, as complemented by reported d etections of possible 
companion cand idates for SN e 1993J (iMaund et al.ll2004 
iFox et aLll20l4 and 201 Idh (|Folatelli et al.ll2ni4lJ l. 

5. CONCLUDING REMARKS 


In this paper, we have reported the late-time evolu¬ 
tion of Type Hb Supernova (SN Hb) 2013df, for which 
the progenitor candidate had been identified as a yel¬ 
low supergiant. SN 2013df showed a dramatic change in 
its spectra at ^ 1 year. On day 626, it has completely 
lost resemblance to other SNe Hb/Ib/Ic, but has been 
dominated by broad and strong Ha and He I emission 
lines arising from a thin shell expanding at the velocity 
of ~ 10, 000 km s“^. Indeed, the only previous example 
showing this characteristics was SN 1993J. Namely, SN 
2013df turns out to be a twin of SN 1993J, not only in 
the early-phase properties but also in the late-time prop¬ 
erties. This behavior is quite different from that seen in 
another well-observed SN Hb 201 Idh, which also showed 
noticeable differences in the (very) early-phase properties 
and in the progenitor than SNe 2013df and 1993J. 

We have shown that this late-time evolution is fully 
explained by a change in the energy input from the ®®Co 
decay to the interaction between the SN ejecta and dense 
CSM. We have derived some characteristic quantities in 
the SN ejecta and CSM under the interaction scenario, 
including the steep density gradient in the outermost SN 
ejecta, the mass and composition within the H-rich en¬ 
velope as are consistent with models for SN 1993J (but 
the mass substantially larger than SN 201 Idh), the CSM 
density, and the mass loss in the final centuries toward 
the explosion. 

SN 2013df on day 626 was found to be brighter than a 
canonical Co-powered model prediction at least by two 
magnitudes, as consistent with the SN-CSM interaction 
scenario. A pros is that we have a chance to observe 
the SN-CSM interaction for a long span to investigate 
the mass loss history further toward the past, but a con 
is that it will raise a hurdle for the confirmation of the 
progenitor disappearance and the detection of a possi¬ 
ble companion. For the SN-CSM interaction study, we 
expect that the spectral features will further change so 
that [O HI] 4959, 500 7 will increase in stren gth as ob¬ 
served for SN 1993J (|Matheson et al.l l2000all and also 
for middle-aged SNe Hn and Hp at > 10 years after the 
explosion (jMilisavlievic et ^l2012tl . This transition will 
take place when the RS reaches to the bo ttom of the H- 
rich layer (e.g., iMilisavlievic et aI1l2012fl . and therefore 
observing this transition will tell us further on the simi¬ 
larity and difference between SNe 2013df and 1993J. On 
the progenitor, we predict that, taking into account the 
interaction power, the SN will keep brighter in the V and 
i?-bands than the progenitor candidate until early 2016, 
or even so until late 2017. A continuous follow-up of this 
SN is therefore very interesting in many aspects. 

The mass loss rate corresponding to the derived CSM 
density is (5.4 ± 3.2) x 10~^Mq yr“^ (for the wind ve¬ 
locity of 20 km s“^), which is similar to SN 1993J but 
larger than SN Hb 2011dh by more than an order of 
magnitude. The striking similarity between SNe 2013df 
and 1993J in the progenitor (candidate) and the CSM 
environment, and the contrast in these natures as com¬ 
pared to SN 201 Idh, infer that these is a link between 
the natures of the progenitor and the mass loss, where 
the SNe Hb with a more extended progenitor have ex¬ 
perienced a much stronger mass loss than those with a 
more compact progenitor in the final centuries toward the 
explosion. We note that a similar conclusion has been in- 





















dependently reac hed by a detailed mo deling; of radio and 
X-ray properties (iKamble et al.l[2?)T^ . strengthening the 
case. Also, there is accumulating evidence that a binary 
interaction is a key in the whole SN Ilb population in¬ 
cluding the present study, and so far there is no single SN 
lib for which a single sta r evolution is favored (see, e.g., 
iKuncaravakti et ^1201511 . The link suggests that there 
might be two classes in SNe lib in terms of characteris¬ 
tic binary evolution schemes: it might indicate that the 
timing between the strong binary interaction and the ex¬ 
plosion might be different between the two classes, which 
might be a key to understand a unified scheme of the yet- 
unresolved evolution toward SNe Ilb/Ib/Ic in general. 
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Fig. 1.— V and i2-band images of SN 2013df at 199 and 626 days after the explosion. 





10 


10.0 


(/> 

N 

E 

o 

Q) 

O 


1.0 


0.1 


K 

3 

U- 0.0 


4000 5000 6000 7000 8000 9000 10000 

Rest Wavelength (A) 



Fig. 2.— Time sequence of late-time spectra of SN 2013df from day 155 to day 626 after the explosion. For demonstration, the fluxes 
in the spectra except for the one on day 199 are multiplied by the amou nts indicated in the figure. The extinction is corrected for with 
E{B — V) = 0.1 mag UMorales-Garoffolo et al'1l2014l : [Van Dvk et al.ll2014l‘) . assuming Ry = 3.1. 
















11 



Rest Wavelength (A) Wavelength (A) 


Fig. 3.— (a) Late -time spectra of SN 2013df (red) as compared to those of SNe 1993J (black: IMatheso n et al.]|2000a|[ hl) and 2 011dh 
('^rav: [ S hivvers et al.l f2Q13ll at the similar epoch s. The exti nction is corrected for with E(B — V) = 0.1 ([Morale^^G^offolo et al.ll2Ql^ : 
IVan Dvk et al.ll2014li . 0.17 l|Ergon et al.l 1201^ 1. and 0.07 l|Ergon et al.l 12014^ mag for SNe 2013df, 1993J, and 2011dh, respectively, 
assuming Ry = 3.1. The flux scale on the vertical axis is for SN 2013df at day 119; the flux at day 199 is added by a constant, and that 
at day 626 is multiplied by a factor of 15 for presentation. The fluxes of the comparison spectra are either multiplied or added by artificial 
amounts for presentation. The spectrum of SN 2013df at day 626 is smoothed with a 5 pixel boxcar filter, (b) A close comparison between 
the spectra of SN 2013df at day 636 (non-smoothed spectrum shown by gray, and the smoothed one shown by red) and that of SN 1993J 
at day 670 (black). 
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Fig. 4.— (a) Ha and (b) He I 5876 (or Na ID) in the spectrum of SN 2013df at day 626 (red) as compared to those of SN 1993J (black). 
In (a), we plot expected profiles computed for the emission from a thin shell with the outer edge velocity of I4h = 10,000 km s~^ and 
the shell width {^Rsh/^sh) 0.2 (green-dashed) and 0.3 (magenta-dashed): The thin shell predicts a flat profile below the wavelength 
corresponding to I4h ~ while a parabola between f/gh ~ 14;^, where AVgh = ARsYi/t. The former (Ai?sh/-Rsh = 0.2) provides 

a reasonable fit to the observed profile, while the latter (Ai?sh/-Rsh = 0.3) is clearly rejected. In (b), Ha of SN 2013df at day 626, as 
divided by a factor of 3.2, is shown (blue). The rest wavelengths of possible contaminating lines are denoted by dashed vertical lines ([O I] 
5577, [N II] 5754, [N H] 6548, [N II] 6583). (c) Evolution of SN 2013df from day 155 to 626, where the spectrum on day 626 multiplied by 
a factor of 20, 8, and 2.5 is compared to the spectra on day 155, 199, and 333, respectively, (d) An expanded view on the evolution of the 
Ha profile from day 199 (gray) to day 626 (red). The thin shell models that provide a reasonable fit to each of the spectra are shown for 
the one at day 199 (green) and day 626 (magenta), which quantify the velocity decrease between the two epochs. 
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Fig. 5. — V and /?-band light curves of SN 2013d f (red points: FMorales-GarofFo lo e^ay[2^^ ^nd this paper) as compared to those of 
SNe 1993J (black open squares: IBarbon et al.iri995f) and 2011dh fgrav filled squares: lErgon et al.ll2014bfl . The magnitudes of SNe 1993J 
and 2011dh are artificially shifted to match to the magnitude scale of SN 2013df. The black-dashed line shows the ^®Co decay line (for a 
full deposition). In the /?-band, the contribution of the interaction power and its evolution (stars) are estimated by the Ho; flux for epochs 
where the spectra are available (Figure 4c). The same behavior is assumed for the V-band. A fit to the interaction light curves is shown for 
a functional form of exponential (magenta-solid) and power-law (magenta-dashed). The ^®Co decay power (represented by the late-time 
evolution of SN 2011dh) is insufficient to explain the brightness of SN 2013df on day 626. Also, the late-time LC evolution of SN 2013df is 
much flatter than the prediction from the ^®Co-decay input, requiring another power source. The inserted panels show the ^®Co decay (full 
deposition) and extrapolation of the interaction light curves in the later epochs, together with the magnitude of the candidate progenitor. 
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Fig. 6.— The relations (a) between the pseudo-bolometri c optical-NIR lum inosity in the early ‘cooling’ phase (upper limits for SN 
Ilb 2008ax and for SN Ib iPTFlSbvn: IPastorello et al1l2008l ; ICao et ^120131) and the CSM density, (b) between the progenitor radius 
and the CSM density, and (c) between the progenitor radius and the mass loss rate. The CSM density is described by A* as defined by 
PCSM = 5 X To convert the CSM density to the mass loss rate, we have assumed the wind velocity of 20 km s~^ for SNe lib 

and 1000 km s“^ for SNe Ib/c. The solid estimate of the CSM density has been provided for SNe 1993J l|Fransson et al.iri996fl and 2011dh 
JMaeda et al.ll2 014all using the com bined information from radio and X-ray data. For SN 2008ax, we conv ert the value estimat ed with radio 
data JChevalier &: Soderberg|l201Cf) to fit to the non-thermal acceleration e fficiency derived for SN 2011d h JMaeda et al.l2014alh The typical 
CSM densities assumed for SNe Ib/c are in dicated by th e shaded region JChevalier &: FYanssonl [20061) . while it may contain a systematic 
uncertainty of up to an order of magnitude JMaeda||20 l'^. The pr ogenitor radii are plotted in (b, c) for those with the prog enitor detection 
(SNe 1993J and 2011dh: IMaund et aL||2004|; |Van Dyk et aL||2013) and with the progenitor candidate detection (SN 2013df:|Van Dyk et aL| 
|20T3). For SNe Ib/c, upper limits on the magnitudes of some SN Ib/c sites in pre-SN images are consistent with a WR star, including a 
possible progenitor detection of iPTF13bvn, i.e., < lOR© as inferred by the shaded area. 




























